Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1788 1438 (N45W48 L=139)  1270    143        0 (1725/1394-1846/1461)   8.5 W38-W59 N41-N47
      2  546 1206 (S09E33 L= 58)   224     12        1 ( 527/1187- 563/1227) -45.0 E34-E31 S09-S07
      3 1738 1078 (N18W43 L=135)   600    155        1 (1700/1010-1789/1130) -16.7 W39-W49 N16-N21
      4  497  996 (S22E31 L= 60)  1031    412        0 ( 314/ 954- 627/1035)  27.6 E46-E22 S28-S16
      5  296  761 (S41E48 L= 43)  1142    161        1 ( 252/ 697- 348/ 838)  42.5 E56-E41 S46-S35
      6 1298  527 (S25W33 L=125)  1243    231        4 (1210/ 467-1391/ 569)  61.4 W27-W38 S31-S20
      7 1011  343 (S46W27 L=118)  1327    217        2 ( 976/ 286-1089/ 441) -15.5 W18-W36 S50-S40
      8 1178  386 (S37W35 L=126)   515     78        1 (1140/ 357-1213/ 417) -14.0 W30-W38 S38-S35
  

    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