Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  357 1576 (N10E63 L=254)  2478    358        8 ( 307/1464- 402/1698) -62.1 E68-E57 N02-N19
      2  422 1623 (N14E59 L=258)   567    195        3 ( 388/1596- 479/1645)  23.2 E62-E53 N12-N15
      3  418 1527 (N08E54 L=263)   351     70        0 ( 381/1495- 441/1556)  53.1 E55-E51 N05-N10
      4  859 1413 (N11E19 L=297)  1005    340        2 ( 773/1302- 970/1545) -29.1 E28-E10 N05-N16
      5  187 1146 (S18E61 L=255)  1988    559        6 (  85/1081- 330/1203)  20.6 E72-E48 S21-S12
      6  773 1079 (S10E16 L=301)   258    103        2 ( 735/1047- 799/1132) -78.7 E17-E13 S12-S07
      7  364 1011 (S23E42 L=275)  2615    817        3 ( 121/ 929- 614/1121)  13.4 E65-E23 S31-S15
      8  774  903 (S20E11 L=305)   429    105        0 ( 720/ 865- 817/ 931)  -7.1 E15-E07 S21-S18
      9 1045  839 (S17W06 L=323)   222     90        0 ( 995/ 806-1085/ 881) -15.9 W02-W09 S18-S15
     10 1518  712 (S10W38 L=355)   803    219        2 (1455/ 637-1573/ 796)  85.2 W35-W41 S16-S04
     11  569  723 (S36E22 L=295)   305     93        1 ( 525/ 704- 615/ 745)   0.0 E25-E17 S36-S34
     12 1444  353 (S32W52 L=  8)  2748    595        3 (1245/ 302-1634/ 400)  20.6 W35-W68 S36-S24
  

    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