Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  701 1698 (N39E39 L=108)   504     77        0 ( 664/1681- 750/1711)   7.1 E42-E34 N37-N39
      2 1098 1469 (N29E02 L=145)  2231    491        2 ( 897/1360-1294/1563) -15.4 E17-W12 N24-N32
      3 1285 1233 (N16W13 L=160)   267     77        0 (1256/1210-1326/1254)  11.3 W11-W16 N14-N17
      4  308 1184 (N00E51 L= 96)   911    120        1 ( 265/1130- 346/1238)  56.3 E54-E48 S03-N03
      5  630 1183 (N04E26 L=120)   219     34        1 ( 611/1161- 649/1203) -45.0 E27-E25 N03-N05
      6 1856 1118 (N17W61 L=208)   588     64        0 (1828/1095-1891/1136)  12.5 W57-W66 N15-N18
      7 1337  745 (S12W23 L=170)   681    202        0 (1250/ 701-1409/ 807)  41.2 W18-W27 S15-S07
      8  684  756 (S21E18 L=129)   224     39        1 ( 659/ 735- 706/ 781)  45.0 E19-E16 S21-S18
      9 1426  597 (S20W33 L=180)  2148    536        0 (1213/ 530-1581/ 747)  40.4 W20-W40 S27-S08
     10 1241  257 (S48W37 L=184)   396     66        1 (1206/ 244-1276/ 274)   0.0 W33-W41 S49-S46
  

    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