Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1340 1686 (N35W47 L=149)  3358    538        4 (1197/1526-1546/1762)  -5.7 W30-W70 N27-N39
      2  593 1127 (N23E24 L= 79)   330     62        1 ( 571/1097- 619/1162) -90.0 E25-E22 N21-N25
      3 1525  952 (S11W27 L=130)   252     57        1 (1503/ 929-1552/ 980) -63.4 W26-W28 S12-S08
      4 1731  894 (S21W42 L=145)  1085    250        1 (1639/ 857-1803/ 964)  -7.6 W34-W49 S23-S15
      5 1812  774 (S32W50 L=153)  2951    386        3 (1736/ 627-1925/ 906)  22.9 W40-W68 S39-S26
      6  170  803 (N12E71 L= 32)  3030    185        3 ( 112/ 739- 219/ 870) -29.7 E80-E64 N07-N15
      7 1421  753 (S19W16 L=119)   416    128        0 (1403/ 693-1442/ 812)  60.3 W13-W18 S22-S15
      8  373  570 (S05E56 L= 47)   583     91        0 ( 324/ 553- 414/ 591)   0.0 E61-E51 S06-S03
      9 1269  471 (S33E01 L=102)   377     79        1 (1242/ 435-1289/ 517)  59.0 E03-W01 S34-S29
     10  755  351 (S27E38 L= 64)   370     43        0 ( 725/ 327- 783/ 375)  18.4 E41-E35 S27-S25
  

    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