Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  856 1505 (N28E18 L= 52)  1378    377        5 ( 748/1381- 965/1628) -33.7 E29-E08 N21-N35
      2 1177 1598 (N40W03 L= 73)   254     22        1 (1158/1578-1200/1616) -18.4 W01-W04 N38-N41
      3  527 1419 (N18E39 L= 31)  2261    469        8 ( 440/1280- 586/1563) -40.1 E51-E31 N10-N26
      4 1293 1241 (N17W14 L= 85)   900    303        2 (1198/1179-1356/1313)  36.9 W08-W18 N13-N21
      5 1519 1216 (N18W30 L=101)   272     90        0 (1506/1167-1537/1265) -71.6 W29-W31 N15-N21
      6  866 1124 (N05E11 L= 60)   241    106        0 ( 808/1113- 918/1137)  15.9 E14-E07 N03-N05
      7 1763 1070 (N12W51 L=122)   886    178        2 (1728/ 986-1800/1137)  68.2 W48-W54 N06-N16
      8  228  976 (S12E57 L= 14)  1093    238        1 ( 166/ 953- 297/1017)  -7.1 E65-E49 S13-S10
      9 1329  399 (S36W33 L=104)  1772    356        4 (1195/ 371-1483/ 436)  12.8 W22-W46 S39-S31
  

    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