Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  548 1779 (N26E60 L=128)  1117    249        3 ( 524/1731- 571/1820) -28.6 E65-E54 N23-N29
      2  917 1800 (N38E35 L=153)   405     78        0 ( 877/1782- 954/1819)   9.5 E38-E32 N36-N38
      3 1176 1651 (N35E11 L=177)   261     64        1 (1143/1636-1212/1674)   9.5 E13-E07 N33-N35
      4 1741 1036 (N15W44 L=231)   385    102        0 (1715/ 989-1759/1080) -26.6 W40-W46 N13-N17
      5  397 1007 (S23E39 L=149)   290     63        1 ( 374/ 977- 426/1031)  76.0 E40-E36 S25-S21
      6  441  766 (S37E32 L=156)   322     65        0 ( 427/ 732- 457/ 798)  76.0 E32-E30 S38-S34
      7 1368  649 (S17W31 L=219)  1532    378        0 (1210/ 626-1547/ 683)  25.5 W20-W41 S21-S09
      8 1400  509 (S23W40 L=228)   934    262        1 (1321/ 442-1453/ 549)  10.3 W33-W44 S28-S20
      9 1070  266 (S49W36 L=224)   504     51        2 (1038/ 245-1106/ 284)  36.9 W33-W38 S51-S47
     10 1275  236 (S43W51 L=239)   398     67        0 (1248/ 218-1294/ 261)  90.0 W48-W52 S44-S41
  

    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