Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  568 1132 (N23E26 L=318)   545    200        0 ( 481/1105- 658/1160) -24.8 E32-E19 N20-N26
      2 1568  996 (S08W32 L= 16)   446    138        0 (1530/ 938-1622/1024)  14.0 W28-W36 S10-S06
      3  878  866 (N02E12 L=332)   725    188        2 ( 858/ 759- 901/ 959)  65.6 E15-E09 S04-N07
      4 1408  843 (S13W18 L=  3)   803    208        0 (1348/ 750-1445/ 944)  45.0 W12-W21 S16-S07
      5 1707  826 (S23W39 L= 24)   846    230        0 (1655/ 715-1731/ 943)  58.4 W33-W42 S28-S15
      6 1264  721 (S16W07 L=351)   306     91        0 (1223/ 697-1315/ 743)   0.0 W03-W10 S17-S14
      7 1378  596 (S27W11 L=356)  1201    353        0 (1215/ 473-1504/ 731)  18.4 E02-W22 S31-S21
      8 1734  669 (S34W41 L= 25)   599     81        0 (1693/ 640-1770/ 701)   7.1 W36-W44 S35-S32
      9 1038  495 (S25E13 L=332)   266     61        1 (1010/ 466-1071/ 525)  18.4 E15-E09 S25-S23
     10  976  275 (S39E26 L=318)   476    125        1 ( 954/ 233-1002/ 319)  36.9 E30-E22 S42-S36
  

    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