Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  909 1693 (N38E26 L=159)   486     93        0 ( 871/1667- 956/1719)  -7.1 E29-E21 N36-N39
      2  473 1618 (N23E54 L=131)   374     96        0 ( 456/1590- 493/1648) -29.7 E57-E50 N21-N25
      3 1235 1527 (N34W03 L=188)  2046    448        5 (1043/1427-1472/1631)  -5.0 E13-W21 N30-N37
      4  496 1496 (N16E46 L=139)   350     92        0 ( 480/1453- 513/1540) -71.6 E47-E43 N13-N19
      5  759 1417 (N16E25 L=160)   245     64        0 ( 744/1380- 781/1454)  90.0 E25-E23 N13-N18
      6  527 1216 (S00E35 L=150)   621    125        1 ( 475/1174- 582/1263)  50.2 E37-E31 S03-N03
      7  951  894 (S11E02 L=183)   219     58        0 ( 926/ 866- 985/ 924) -33.7 E03-W00 S11-S09
      8 1115  858 (S10W09 L=194)   840    216        4 (1023/ 822-1224/ 907)  32.5 W03-W14 S12-S04
      9 1460  563 (S20W39 L=224)   597    121        1 (1396/ 539-1518/ 589)  35.5 W35-W42 S22-S17
     10 1108  523 (S30W17 L=202)   577    123        4 (1047/ 498-1175/ 558)   6.3 W12-W21 S31-S28
  

    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