Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1287 1425 (N21W27 L= 98)   413    115        2 (1257/1371-1321/1483)  29.7 W24-W31 N17-N23
      2 1843 1278 (S05W64 L=135)  2633    241        3 (1795/1186-1874/1354)  68.2 W58-W70 S10-N00
      3 1126 1212 (N14W11 L= 82)   281     97        0 (1070/1193-1177/1244)  -8.1 W07-W14 N12-N15
      4  655 1086 (N19E20 L= 51)  1020    304        2 ( 530/1048- 809/1142) -34.5 E27-E10 N13-N24
      5 1631 1048 (S10W36 L=107)   420     86        0 (1615/ 990-1650/1095)  51.3 W34-W38 S12-S07
      6  189  865 (N16E63 L=  9)  1156    227        6 ( 159/ 796- 205/ 954)  48.4 E67-E58 N12-N21
      7 1612  783 (S24W29 L=100)   967    192        1 (1573/ 704-1637/ 881)  59.0 W25-W32 S28-S18
      8  670  711 (S02E29 L= 42)   223     50        0 ( 654/ 680- 684/ 739) -90.0 E29-E27 S04-N00
      9 1365  585 (S29W08 L= 79)  1194    291        1 (1249/ 518-1435/ 698)  11.3 W00-W15 S32-S23
     10  834  632 (S11E22 L= 50)   376    146        0 ( 764/ 600- 902/ 669)   0.0 E26-E16 S11-S09
     11  910  517 (S20E21 L= 50)   547    196        0 ( 888/ 441- 938/ 615)  45.0 E25-E16 S23-S14
  

    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