Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1265 1725 (N39W40 L=201)  1049    207        1 (1167/1675-1361/1760)  -2.9 W30-W50 N37-N40
      2 1428 1535 (N23W40 L=201)   319    101        0 (1390/1508-1473/1566) -63.4 W37-W41 N19-N25
      3  481 1336 (N34E30 L=132)   770    192        0 ( 396/1316- 577/1354) -10.6 E37-E21 N31-N35
      4 1107 1266 (N15W10 L=171)   276     89        0 (1059/1246-1148/1281) -31.0 W07-W12 N12-N16
      5 1808 1213 (S05W57 L=218)   427     84        0 (1782/1167-1827/1249) -78.7 W55-W58 S07-S02
      6  680 1227 (N22E16 L=145)   358     80        0 ( 640/1199- 724/1248)   9.5 E19-E13 N21-N23
      7 1739 1048 (S12W46 L=207)  1292    273        2 (1672/ 971-1796/1130) -61.4 W41-W52 S17-S06
      8  278 1008 (N18E51 L=111)  1650    326        1 ( 145/ 943- 388/1108)   6.1 E67-E39 N14-N21
      9 1839 1026 (S16W56 L=217)   470     65        0 (1820/ 995-1858/1070)  45.0 W53-W59 S18-S14
     10 1433  947 (S10W22 L=183)   223     53        0 (1414/ 922-1448/ 980)  45.0 W20-W23 S11-S08
     11 1924  913 (S26W68 L=229)   631     50        1 (1903/ 885-1949/ 935)   6.3 W64-W74 S27-S23
     12 1560  798 (S22W28 L=189)   347     61        0 (1539/ 766-1577/ 837)  68.2 W26-W29 S24-S19
     13 1388  729 (S22W15 L=176)   773    204        0 (1322/ 647-1440/ 817)  32.5 W09-W20 S25-S18
     14 1184  664 (S21W01 L=162)  1013    287        0 (1146/ 546-1252/ 804)  52.4 E03-W08 S27-S14
     15  870  651 (S15E18 L=143)   301    105        0 ( 840/ 614- 912/ 699)  29.7 E20-E13 S16-S12
     16 1412  597 (S31W15 L=176)   515    118        0 (1356/ 568-1460/ 639)   7.1 W10-W18 S32-S29
     17  358  582 (S08E56 L=105)   343     32        1 ( 341/ 561- 378/ 604) -45.0 E58-E54 S09-S07
     18  824  446 (S26E28 L=133)   553    151        1 ( 787/ 376- 850/ 517)  38.7 E33-E23 S30-S22
     19 1630  449 (S50W35 L=196)   957    117        0 (1565/ 426-1687/ 479)  -3.6 W27-W43 S52-S48
  

    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