Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  460 1587 (N28E46 L=273)   570    183        1 ( 394/1569- 517/1610)  -9.5 E52-E40 N26-N29
      2 1587 1463 (N27W36 L=355)   527    150        1 (1513/1445-1646/1477)   5.7 W30-W40 N26-N28
      3 1316 1297 (N14W16 L=335)   216     60        1 (1303/1264-1331/1341)  90.0 W15-W16 N11-N16
      4  124  893 (S16E69 L=250)  3297    495        9 (  85/ 723- 178/1009) -80.0 E75-E61 S25-S08
      5  291  780 (S23E51 L=268)   998    271        1 ( 241/ 677- 322/ 890)  72.9 E56-E47 S29-S16
      6  252  803 (S22E54 L=265)   481    106        2 ( 221/ 774- 279/ 844)  45.0 E57-E51 S23-S19
      7  700  788 (S21E19 L=300)   404     76        1 ( 660/ 761- 734/ 825)  36.9 E21-E16 S22-S18
      8 1841  672 (S17W65 L= 24)   833    138        1 (1803/ 639-1877/ 709)  35.5 W60-W69 S19-S14
      9 1638  601 (S25W47 L=  7)   577    106        0 (1592/ 559-1681/ 657)  60.3 W44-W49 S27-S20
     10  562  613 (S33E30 L=289)   273     85        0 ( 540/ 588- 598/ 643) -14.0 E31-E27 S34-S30
     11 1000  538 (S35W02 L=321)   251    142        0 ( 938/ 482-1063/ 572)  45.0 E01-W06 S39-S32
     12 1638  446 (S34W54 L= 13)  3351    547        6 (1459/ 377-1813/ 515)  30.0 W41-W68 S41-S26
  

    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