Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1583 1528 (N26W45 L=  1)   866    153        0 (1501/1488-1645/1555) -20.0 W38-W49 N22-N27
      2  401 1201 (N13E40 L=277)   282     72        0 ( 366/1188- 442/1212)   0.0 E42-E36 N12-N13
      3 1450 1031 (S04W26 L=343)   252     62        0 (1423/1009-1488/1049)   0.0 W24-W28 S05-S02
      4 1797  960 (S10W54 L= 10)   769    175        0 (1753/ 894-1846/1042)  38.7 W49-W59 S13-S05
      5  918  806 (S14E08 L=308)   470    194        0 ( 853/ 729- 966/ 887) -61.4 E11-E05 S19-S08
      6  401  769 (S12E44 L=273)   277     60        1 ( 388/ 738- 415/ 803)  76.0 E44-E42 S14-S10
      7  429  662 (S19E44 L=272)   929    272        3 ( 400/ 557- 463/ 771)  52.4 E49-E38 S25-S12
      8 1580  622 (S30W38 L=354)   336     71        0 (1546/ 598-1612/ 647)  31.0 W34-W40 S31-S28
      9  773  529 (S31E22 L=295)   303     55        0 ( 746/ 505- 801/ 552)  26.6 E23-E19 S32-S29
     10  743  377 (S41E29 L=287)   457     52        1 ( 720/ 346- 766/ 402)  45.0 E31-E26 S43-S39
     11  676  276 (S49E43 L=274)   516     79        1 ( 645/ 246- 718/ 300)  13.0 E49-E36 S50-S46
  

    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