Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1144 1884 (N55E31 L= 42)   659     42        1 (1123/1859-1175/1905)   8.1 E35-E27 N51-N57
      2  571 1672 (N19E51 L= 21)  1119    134        1 ( 540/1615- 605/1733) -28.3 E58-E45 N15-N22
      3  207 1171 (S18E62 L= 10)  1184    274        1 ( 149/1080- 285/1287)  51.8 E69-E56 S23-S09
      4 1209 1235 (N10W04 L= 77)   261     81        2 (1189/1188-1227/1272) -59.0 W03-W06 N07-N12
      5 1302 1198 (N11W10 L= 83)   264    114        1 (1286/1140-1330/1240) -68.2 W08-W11 N07-N13
      6  313 1171 (S16E50 L= 22)   425    124        1 ( 297/1111- 326/1231) -60.3 E53-E48 S19-S12
      7 1598 1122 (N16W31 L=103)   236     73        0 (1585/1084-1613/1159) -63.4 W29-W32 N13-N17
      8  199 1044 (S25E61 L= 12)  1875    180        5 ( 154/ 967- 244/1129)  59.0 E67-E55 S30-S20
      9  386 1028 (S23E41 L= 32)   972    286        0 ( 269/ 968- 494/1091)  33.7 E51-E33 S28-S16
     10 1064 1002 (S07W03 L= 75)   317     82        1 (1034/ 968-1109/1030)  11.3 W00-W05 S09-S05
  

    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