Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1502 1631 (N25W54 L= 38)   784    101        1 (1472/1591-1531/1674)  56.3 W52-W57 N22-N28
      2 1500 1506 (N19W45 L= 29)   421     87        0 (1474/1466-1524/1548)  20.6 W41-W49 N17-N20
      3 1267 1486 (N25W29 L= 13)   394    121        0 (1205/1464-1319/1511) -50.2 W26-W32 N22-N28
      4 1141 1378 (N23W17 L=  1)   759    280        2 (1020/1329-1244/1416) -30.3 W08-W21 N17-N26
      5  933 1302 (N25W02 L=346)   221    101        0 ( 912/1267- 962/1351)  31.0 W00-W05 N22-N26
      6  869 1282 (N25E02 L=341)   223     80        1 ( 840/1251- 888/1308) -76.0 E03-W00 N23-N27
      7  342 1267 (N38E44 L=300)   473     45        1 ( 315/1240- 363/1290) -45.0 E46-E41 N35-N39
      8  419 1228 (N34E36 L=308)   493     79        2 ( 383/1201- 461/1250) -23.2 E39-E32 N31-N35
      9  893 1204 (N20E03 L=341)   302    121        0 ( 844/1175- 943/1236) -51.3 E05-W00 N17-N23
     10  530 1133 (N25E28 L=315)   410    130        1 ( 500/1089- 572/1190) -74.1 E29-E25 N21-N28
     11  159 1020 (N26E68 L=276)   676     80        0 ( 144/ 983- 175/1059) -76.0 E70-E64 N23-N27
     12  185  933 (N20E64 L=279)  2343    178        6 ( 135/ 847- 227/1035)  61.4 E72-E57 N14-N25
     13 1420  862 (S13W18 L=  2)   989    323        0 (1357/ 739-1482/1008)  40.6 W11-W25 S18-S06
     14  660  897 (N09E24 L=319)   520    164        0 ( 602/ 836- 721/ 959) -68.2 E26-E21 N03-N13
     15 1706  878 (S22W39 L= 23)   691    107        1 (1685/ 819-1725/ 935)  63.4 W37-W41 S25-S18
     16  371  876 (N14E46 L=298)   313     94        1 ( 322/ 853- 401/ 894) -26.6 E49-E43 N11-N14
     17 1174  612 (S21E03 L=341)   345    151        0 (1137/ 587-1232/ 656) -45.0 E04-W01 S23-S17
     18 1416  589 (S30W12 L=355)   629    123        1 (1374/ 526-1458/ 631)  32.0 W07-W15 S33-S26
     19  927  258 (S38E34 L=309)   601    120        1 ( 899/ 202- 950/ 317)  57.5 E39-E31 S43-S32
  

    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