Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  720 1845 (N36E53 L=169)   629    177        0 ( 701/1815- 741/1880) -24.0 E57-E48 N33-N38
      2  609 1291 (S02E31 L=190)   219     68        1 ( 582/1275- 642/1312)   0.0 E33-E29 S03-N00
      3  363 1154 (S15E45 L=176)   532    113        0 ( 319/1109- 410/1189)  50.2 E47-E41 S18-S12
      4 1850 1153 (N27W52 L=273)   491     79        1 (1827/1123-1870/1184) -36.9 W49-W54 N25-N29
      5 1696 1089 (N18W39 L=260)   299     65        0 (1673/1052-1713/1122) -36.9 W36-W41 N15-N19
      6 1574  989 (N08W33 L=254)   882    179        0 (1518/ 912-1636/1068) -24.4 W27-W38 N05-N11
      7  482  998 (S22E32 L=189)   311    103        0 ( 453/ 966- 526/1048) -26.6 E35-E28 S22-S19
      8 1890  891 (N14W64 L=285)  3578    417        2 (1852/ 742-1922/1025) -43.0 W54-W73 N06-N21
      9 1058  937 (S10W04 L=226)   374    153        0 (1045/ 855-1069/1013) -53.1 W01-W07 S14-S06
     10  693  920 (S21E16 L=205)   212     72        0 ( 656/ 907- 733/ 935)  11.3 E18-E13 S21-S19
     11  343  742 (S41E40 L=181)  1055    326        3 ( 247/ 716- 460/ 766)  15.9 E50-E29 S44-S36
     12  949  630 (S30W08 L=229)   562    287        0 ( 802/ 598-1083/ 692)  12.5 E03-W15 S31-S26
     13  802  430 (S48W06 L=227)   265     46        0 ( 780/ 411- 827/ 449) -11.3 W03-W08 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