Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  662 1619 (N38E35 L=171)  1300    279        2 ( 617/1477- 715/1722) -42.0 E46-E26 N28-N46
      2 1082 1483 (N32W01 L=208)  1405    447        1 ( 954/1330-1233/1570) -34.4 E09-W11 N22-N37
      3  412 1374 (N19E47 L=160)   396     88        0 ( 394/1338- 443/1418) -29.7 E50-E43 N17-N22
      4 1779 1164 (N15W55 L=262)  1034    189        0 (1751/1094-1807/1246) -90.0 W52-W57 N10-N19
      5 1850 1084 (N10W63 L=270)   452     67        1 (1835/1050-1867/1119)  78.7 W60-W65 N07-N12
      6 1754  857 (S04W53 L=260)   484     32        1 (1730/ 833-1778/ 879) -11.3 W50-W55 S05-S02
      7  978  756 (S15E01 L=206)   230     73        2 ( 960/ 719- 994/ 791)  90.0 E02-W00 S17-S12
      8 1855  702 (S14W72 L=279)  2480    175        5 (1826/ 629-1889/ 768)  77.5 W67-W80 S18-S09
      9 1756  649 (S18W60 L=267)   683     43        2 (1730/ 629-1784/ 671) -14.0 W56-W64 S18-S16
     10 1229  577 (S25W17 L=224)   283     62        0 (1193/ 564-1266/ 594)  21.8 W14-W19 S25-S23
     11  906  493 (S33E05 L=202)   275     64        0 ( 871/ 478- 942/ 510)  21.8 E07-E01 S34-S31
     12 1091  425 (S37W10 L=217)   328    114        0 (1039/ 407-1148/ 444)  20.6 W06-W14 S38-S35
  

    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