Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  570 1676 (N29E52 L= 70)   429     85        0 ( 548/1654- 603/1714) -18.4 E56-E47 N28-N32
      2  829 1422 (N18E21 L=102)   783    209        3 ( 769/1347- 886/1502) -41.6 E26-E15 N14-N22
      3  495 1166 (S03E36 L= 87)   542    117        0 ( 440/1146- 537/1193)  26.6 E39-E32 S05-S01
      4 1861 1011 (N15W62 L=185)  1247    235        3 (1820/ 930-1894/1107) -51.3 W56-W68 N10-N20
      5  927  955 (S07E04 L=118)   598    243        1 ( 853/ 878-1005/1039)  63.4 E07-E01 S13-S01
      6 1240  828 (S09W17 L=139)   233     96        0 (1222/ 790-1264/ 869)  80.5 W15-W17 S11-S05
      7 1250  752 (S13W19 L=142)   299    115        1 (1225/ 704-1284/ 803)  90.0 W18-W21 S16-S09
      8 1158  700 (S18W15 L=137)   246    120        0 (1102/ 666-1205/ 726)   0.0 W11-W18 S19-S16
  

    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