Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1605 1650 (N35W58 L=296)   544    105        0 (1572/1634-1641/1664) -26.6 W55-W61 N33-N36
      2 1222 1371 (N24W18 L=255)   372     94        0 (1180/1328-1267/1413) -50.2 W15-W20 N21-N27
      3  751 1184 (N19E15 L=222)  1587    396        5 ( 561/1095- 922/1277) -31.3 E28-E05 N11-N26
      4  568 1028 (N11E29 L=208)   607    171        0 ( 509/ 948- 629/1073) -48.8 E33-E25 N05-N13
      5  238  971 (N09E58 L=180)   358     78        1 ( 218/ 934- 253/ 998)  63.4 E59-E55 N07-N11
      6 1784  802 (S20W52 L=289)  2702    317        4 (1683/ 666-1847/ 936)  41.2 W42-W59 S27-S12
      7 1164  821 (S08W06 L=244)   260     25        0 (1142/ 794-1188/ 846)  45.0 W04-W07 S09-S06
      8 1233  745 (S14W09 L=247)   860    334        0 (1086/ 651-1333/ 836)  25.0 W00-W16 S17-S08
      9 1053  775 (S10E01 L=236)   191     95        0 (1024/ 736-1064/ 817)  59.0 E03-W00 S11-S06
     10 1579  657 (S25W34 L=271)   326     98        1 (1549/ 607-1595/ 699)  71.6 W31-W34 S28-S22
     11  278  615 (S13E62 L=175)   718    105        1 ( 264/ 571- 301/ 662)  63.4 E64-E59 S16-S10
     12 1478  614 (S26W26 L=263)   363     72        1 (1447/ 569-1507/ 646)  45.0 W23-W27 S29-S24
  

    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