Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1263 1623 (N35E01 L= 34)   319     59        1 (1240/1595-1292/1651) -26.6 E03-W01 N33-N37
      2 1364 1325 (N18W12 L= 47)   235     77        0 (1330/1288-1393/1360) -21.8 W09-W14 N16-N19
      3  197 1146 (S17E60 L=335)   992    207        1 ( 171/1049- 223/1256) -85.2 E62-E58 S22-S10
      4  426 1146 (S13E40 L=355)  1025    246        2 ( 390/1053- 464/1243) -85.2 E43-E36 S19-S07
      5  492 1163 (S11E35 L=360)   294     49        0 ( 470/1132- 519/1193)  78.7 E36-E33 S13-S08
      6 1890  987 (N16W60 L= 95)  2560    344        4 (1835/ 919-1932/1061) -53.1 W54-W67 N11-N20
      7  868  935 (S16E07 L= 28)   454    161        0 ( 815/ 864- 905/1008) -29.1 E11-E02 S19-S12
      8 1516  879 (S02W32 L= 67)  1052    217        2 (1424/ 831-1611/ 959)  -3.8 W24-W39 S04-N00
      9 1033  819 (S18W06 L= 41)   212     28        0 (1018/ 792-1048/ 846) -71.6 W04-W06 S19-S16
     10  537  788 (S32E26 L= 10)  1043    322        0 ( 403/ 761- 685/ 815)  12.8 E36-E14 S35-S29
     11 1864  771 (N03W66 L=101)   503     57        0 (1848/ 736-1881/ 802) -76.0 W64-W68 N01-N05
     12 1693  717 (S06W51 L= 86)   496    152        0 (1641/ 678-1728/ 754)   0.0 W45-W55 S06-S04
     13 1591  598 (S15W48 L= 83)  1068    223        1 (1487/ 560-1682/ 650)  34.7 W40-W53 S18-S09
     14 1356  630 (S20W31 L= 66)   235     69        0 (1321/ 610-1384/ 646)   0.0 W28-W33 S20-S19
     15 1738  512 (S14W66 L=101)  3994    436        6 (1676/ 437-1801/ 598) -73.3 W59-W75 S19-S09
     16 1645  487 (S19W58 L= 93)   441    122        1 (1622/ 461-1680/ 515)  -8.1 W55-W62 S19-S17
     17 1270  479 (S31W32 L= 67)   831    136        0 (1195/ 457-1346/ 509)  40.6 W27-W36 S34-S27
     18 1535  465 (S23W51 L= 86)   461    118        3 (1512/ 425-1559/ 504) -68.2 W47-W52 S25-S20
  

    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