Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  504 1534 (N22E41 L=212)  1499    641        1 ( 354/1446- 694/1654) -20.4 E60-E25 N17-N30
      2  400 1135 (S03E41 L=212)   546    219        1 ( 316/1081- 478/1169) -12.1 E48-E34 S06-S01
      3  205  981 (S13E58 L=195)   538    130        2 ( 168/ 940- 237/1012) -29.7 E62-E54 S15-S11
      4  404  944 (S15E40 L=213)   359    105        2 ( 365/ 913- 450/ 976)  33.7 E42-E36 S16-S12
      5  975  856 (S16E01 L=252)  1117    321        4 ( 880/ 742-1122/ 959)  51.3 E06-W06 S23-S08
      6  535  920 (S15E30 L=223)   209     64        1 ( 499/ 897- 570/ 938)  11.3 E32-E27 S16-S14
      7  674  860 (S18E20 L=232)   260     45        2 ( 655/ 835- 697/ 892) -56.3 E21-E18 S19-S16
      8  312  769 (S26E48 L=204)   353    108        1 ( 292/ 719- 341/ 802)  51.3 E50-E45 S28-S23
      9  463  727 (S28E36 L=217)   448    159        0 ( 419/ 670- 513/ 773)  37.9 E40-E31 S31-S24
     10 1245  724 (S21W17 L=270)   285     80        0 (1215/ 686-1266/ 772)  71.6 W15-W18 S23-S17
     11 1618  599 (S24W47 L=300)  1427    214        4 (1532/ 529-1694/ 668)  48.0 W42-W51 S29-S18
  

    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