Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  792 1713 (N33E38 L=269)  1174    201        2 ( 728/1653- 870/1785) -14.7 E48-E29 N30-N37
      2  573 1627 (N22E48 L=258)   828    187        0 ( 525/1556- 601/1688) -33.7 E56-E43 N17-N25
      3 1399 1632 (N46W10 L=316)   691    125        0 (1353/1595-1455/1657) -11.3 W05-W15 N43-N48
      4  628 1544 (N18E40 L=267)   518    109        0 ( 588/1500- 659/1592)  90.0 E40-E38 N14-N21
      5 1777 1191 (N26W49 L=355)  1055    131        3 (1746/1119-1803/1252) -54.5 W45-W52 N22-N29
      6 1321 1176 (N13W14 L=320)   227     76        0 (1299/1142-1346/1214) -36.9 W11-W15 N11-N14
      7 1863 1147 (N26W60 L=  7)   623     98        0 (1832/1106-1897/1189) -11.3 W55-W65 N24-N28
      8 1557  574 (S15W47 L=354)   462    133        3 (1537/ 516-1575/ 643) -41.2 W42-W50 S18-S11
      9 1286  573 (S22W29 L=335)   294    105        0 (1272/ 531-1310/ 626) -39.8 W25-W32 S24-S19
     10 1363  483 (S26W38 L=345)   285     97        0 (1321/ 457-1404/ 510)  63.4 W36-W40 S28-S22
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y