Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1169 1814 (N57E13 L= 36)  1830    240        2 (1063/1768-1266/1850)  43.0 E21-E04 N48-N64
      2 1798 1106 (N20W53 L=102)  2152    263        5 (1733/1024-1858/1186)  84.3 W46-W61 N14-N24
      3  250 1120 (S11E56 L=353)   834    191        1 ( 202/1052- 310/1160)  48.8 E61-E50 S15-S07
      4 1690 1097 (N16W42 L= 92)   931    175        2 (1643/1014-1732/1160) -41.2 W38-W47 N12-N19
      5 1607 1091 (N14W36 L= 85)   242     72        3 (1586/1052-1620/1122) -36.9 W33-W37 N12-N16
      6 1905 1005 (N16W69 L=119)  1289     91        2 (1885/ 955-1935/1053)  90.0 W66-W76 N12-N18
      7 1254  614 (S21W24 L= 73)   278    117        0 (1210/ 579-1300/ 648)  56.3 W21-W25 S24-S18
      8 1720  508 (S17W68 L=117)   584     73        1 (1695/ 476-1752/ 536)  53.1 W65-W70 S19-S15
      9 1499  436 (S27W50 L= 99)  1144    152        3 (1439/ 407-1571/ 475)  40.6 W44-W55 S29-S22
     10 1303  372 (S36W39 L= 88)   487     81        2 (1266/ 354-1349/ 396)   7.1 W34-W42 S36-S33
     11 1169  230 (S50W41 L= 90)   436     48        1 (1145/ 202-1188/ 253) -24.0 W36-W45 S52-S48
  

    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