Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  747 1646 (N27E31 L=147)   747    142        1 ( 718/1592- 796/1694)  90.0 E33-E28 N23-N31
      2  365 1558 (N16E59 L=120)  4812    300       10 ( 288/1432- 433/1650) -38.3 E77-E49 N07-N22
      3  821 1439 (N14E20 L=159)   204     51        0 ( 795/1426- 853/1455)   0.0 E21-E17 N13-N15
      4 1711 1331 (N26W41 L=220)   939    149        3 (1665/1266-1776/1375) -21.8 W37-W47 N23-N28
      5 1894 1212 (N24W63 L=242)  1013    117        3 (1856/1154-1927/1255) -20.0 W57-W68 N20-N26
      6  307 1122 (S11E50 L=129)   421     68        0 ( 280/1096- 340/1153)  36.9 E52-E46 S12-S09
      7  269 1017 (S18E52 L=127)  1601    247        4 ( 215/ 912- 325/1099) -45.0 E58-E46 S23-S12
      8  512  981 (S17E31 L=147)   213     77        0 ( 480/ 951- 543/1011)  45.0 E33-E29 S19-S15
      9 1343  890 (S08W21 L=200)   210     66        0 (1309/ 865-1373/ 913)  45.0 W19-W22 S10-S06
     10  602  750 (S30E23 L=156)   239     33        1 ( 580/ 730- 620/ 778)  56.3 E24-E21 S31-S27
     11  729  735 (S29E14 L=165)   676    193        0 ( 644/ 709- 821/ 773)   0.0 E19-E07 S29-S27
     12  520  740 (S32E29 L=149)   246     64        0 ( 486/ 726- 559/ 759)   0.0 E32-E26 S32-S30
     13 1830  608 (S12W73 L=251)   860     89        3 (1810/ 553-1850/ 643) -50.2 W68-W77 S15-S09
     14  283  600 (S44E59 L=120)   788     80        0 ( 237/ 571- 324/ 631)  12.5 E69-E51 S46-S41
     15 1507  568 (S23W42 L=221)   609    260        0 (1393/ 523-1608/ 626)  39.8 W35-W48 S27-S17
     16  460  433 (S55E42 L=137)  1372    187        3 ( 377/ 385- 559/ 476)  14.0 E63-E27 S59-S49
     17 1618  379 (S31W67 L=245)  3882    281        6 (1484/ 336-1739/ 438)  27.6 W54-W81 S36-S24
     18 1308  294 (S45W46 L=224)   338     42        0 (1282/ 276-1337/ 310)  56.3 W43-W47 S46-S43
  

    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