Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  831 1604 (N36E11 L=184)   381     80        0 ( 784/1584- 873/1626)   9.5 E14-E08 N34-N37
      2 1606 1564 (N30W46 L=241)   453    107        1 (1575/1531-1635/1595)  23.2 W42-W49 N27-N31
      3  766 1303 (N16E15 L=180)   349    118        0 ( 713/1287- 840/1323)  -7.1 E18-E10 N14-N17
      4  151 1237 (N15E65 L=130)  1055    239        0 ( 116/1202- 186/1282) -24.4 E71-E60 N12-N18
      5 1402 1183 (N06W24 L=219)   194     33        1 (1385/1164-1419/1207)  56.3 W22-W24 N04-N07
      6 1072  953 (S07W03 L=198)  1021    296        4 (1036/ 805-1122/1106) -80.0 W00-W05 S15-N02
      7  236  945 (S03E55 L=140)   435    100        0 ( 200/ 914- 270/ 972) -26.6 E59-E52 S04-S01
      8  626  512 (S32E31 L=164)  1020    226        3 ( 518/ 480- 723/ 538)  -3.6 E40-E24 S34-S30
      9  510  449 (S36E43 L=152)   435    102        1 ( 488/ 420- 550/ 478)  15.9 E46-E39 S37-S33
  

    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