Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  561 1541 (N12E43 L=189)  1123    303        1 ( 515/1404- 635/1664) -80.5 E45-E39 N03-N21
      2 1317 1239 (N12W11 L=243)  1975    438        2 (1225/1043-1446/1455) -38.0 W00-W23 N04-N22
      3 1480 1423 (N28W18 L=249)   246     47        0 (1455/1394-1499/1447) -26.6 W15-W19 N26-N29
      4  673 1330 (N02E28 L=204)  1507    391        4 ( 485/1249- 806/1370)   2.5 E41-E18 N00-N04
      5  230 1276 (S09E61 L=171)   431     89        0 ( 201/1248- 254/1310)  45.0 E63-E59 S10-S06
      6  780 1141 (S07E17 L=215)   686    198        0 ( 744/1079- 804/1237) -58.0 E19-E13 S10-S01
      7 1789 1003 (N13W49 L=281)  2839    422        5 (1753/ 799-1820/1213) -59.4 W43-W58 N02-N25
      8  604 1083 (S14E26 L=206)   976    251        0 ( 487/1045- 721/1135)  29.7 E33-E19 S17-S08
      9  975  719 (S25W05 L=237)  1044    320        1 ( 812/ 692-1134/ 764)  35.2 E04-W13 S30-S18
     10  464  713 (S39E31 L=201)  1525    361        4 ( 329/ 696- 606/ 740)  14.0 E43-E19 S42-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