Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1045 1692 (N45W27 L= 75)  1029    153        2 ( 976/1654-1122/1716) -24.4 W20-W33 N43-N49
      2 1396 1560 (N26W41 L= 89)   436     90        2 (1357/1541-1437/1577) -31.0 W38-W44 N23-N27
      3 1590 1455 (N14W49 L= 97)   613    103        0 (1536/1438-1639/1474) -11.3 W44-W54 N12-N15
      4  940 1281 (N23W02 L= 50)   248     62        1 ( 916/1254- 966/1309)  18.4 W00-W03 N22-N25
      5 1844 1146 (S11W59 L=107)  1397    193        2 (1794/1071-1902/1208)  13.2 W51-W69 S14-S07
      6  649 1066 (N18E21 L= 27)   723    210        1 ( 579/1001- 721/1138)  18.4 E27-E15 N16-N21
      7  449 1097 (N25E35 L= 13)   233     45        0 ( 421/1078- 476/1116) -36.9 E37-E32 N23-N26
      8  179  992 (N23E65 L=343)  2125    237        3 ( 155/ 894- 230/1105)  81.3 E68-E58 N17-N30
      9  343  995 (N21E46 L=  2)  1080    252        1 ( 233/ 920- 398/1052) -32.7 E56-E41 N15-N25
     10  988 1021 (N08E02 L= 46)   308     59        2 ( 957/ 991-1018/1051)  33.7 E03-W00 N06-N09
     11  260  892 (N16E56 L=352)   613     94        2 ( 242/ 851- 282/ 949) -90.0 E58-E53 N13-N19
     12 1516  916 (S13W26 L= 74)   262     83        4 (1490/ 883-1548/ 945) -53.1 W24-W27 S14-S10
     13 1760  873 (S24W44 L= 92)   762    151        1 (1720/ 812-1802/ 932)  21.8 W39-W49 S25-S21
     14 1809  814 (S29W49 L= 97)   395     68        0 (1785/ 790-1830/ 851)  21.8 W46-W52 S30-S27
     15 1388  758 (S18W14 L= 62)  1020    344        2 (1226/ 742-1549/ 787) -20.9 W04-W25 S22-S13
     16  517  458 (S14E50 L=358)  1584    407        4 ( 439/ 355- 571/ 636)  30.7 E64-E37 S20-S04
     17  866  467 (S22E26 L= 22)   454    117        2 ( 844/ 409- 890/ 514)  66.8 E27-E23 S25-S18
     18  647  275 (S30E54 L=354)   464     65        0 ( 628/ 238- 669/ 306)  26.6 E59-E48 S32-S27
  

    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