Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1403 1678 (N33W45 L= 75)  1789    302        4 (1270/1621-1524/1728) -38.2 W35-W52 N26-N39
      2  482 1667 (N57E36 L=354)  2135    171        7 ( 423/1595- 542/1723)  73.3 E48-E27 N50-N61
      3 1737 1590 (N20W72 L=102)   785     63        1 (1701/1565-1769/1606)  -3.6 W66-W82 N18-N20
      4  764 1539 (N39E06 L= 23)  1312    257        0 ( 671/1464- 872/1584) -45.0 E13-W00 N31-N43
      5  855 1325 (N23E04 L= 26)   401     84        1 ( 820/1299- 895/1349)  11.3 E06-E01 N20-N24
      6  265 1204 (N27E51 L=339)  1376    268        3 ( 228/1120- 315/1328) -57.3 E57-E45 N21-N35
      7  960 1226 (N15W00 L= 30)   492    116        0 ( 908/1195-1020/1267) -56.3 E02-W03 N11-N18
      8  383 1121 (N20E40 L=350)   880    117        3 ( 345/1060- 429/1174)  63.4 E43-E35 N16-N22
      9 1786 1061 (S11W52 L= 81)  1176    205        0 (1760/ 955-1813/1151)  57.5 W47-W56 S17-S06
     10  303 1020 (N15E48 L=342)  1319    188        4 ( 241/ 947- 358/1077) -37.9 E54-E44 N10-N18
     11  146  999 (N16E67 L=323)  1065    118        2 ( 114/ 957- 192/1051) -33.7 E72-E60 N13-N19
     12 1616  841 (S21W34 L= 64)   602    120        0 (1562/ 800-1685/ 876) -37.9 W30-W39 S24-S17
     13  312  732 (S02E54 L=336)   683    176        0 ( 243/ 681- 365/ 782) -41.2 E59-E49 S05-N02
     14  846  555 (S22E21 L=  9)   406    165        0 ( 803/ 489- 882/ 627)  37.9 E25-E16 S25-S18
     15 1122  564 (S27E04 L= 26)   312     76        1 (1099/ 530-1147/ 601)  68.2 E05-E01 S29-S24
     16  383  478 (S18E60 L=330)  1599    165        8 ( 328/ 419- 421/ 534) -74.1 E65-E53 S22-S14
     17  576  302 (S34E55 L=334)   593     86        2 ( 540/ 277- 612/ 332)   0.0 E60-E50 S35-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