Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1197 1776 (N57W16 L=107)   339     48        1 (1171/1753-1223/1795)  15.9 W12-W19 N55-N59
      2  394 1510 (N31E55 L= 36)   548     86        1 ( 354/1490- 438/1525)   0.0 E60-E49 N29-N31
      3  718 1378 (N24E22 L= 69)   622    202        2 ( 636/1340- 799/1412)  -9.5 E28-E16 N21-N25
      4  870 1251 (N16E11 L= 81)   257    113        0 ( 810/1234- 922/1263)  -8.1 E14-E07 N15-N17
      5 1292 1203 (N15W17 L=108)   337    142        0 (1225/1174-1343/1235) -20.6 W12-W20 N12-N16
      6 1491  748 (S13W32 L=123)   281    174        0 (1428/ 689-1544/ 806)  53.1 W28-W35 S17-S09
      7 1818  731 (S15W63 L=155)   598    100        2 (1787/ 688-1851/ 765)  51.3 W60-W67 S17-S12
      8  736  709 (S18E18 L= 73)   213     99        0 ( 695/ 696- 781/ 728)   9.5 E20-E14 S18-S16
      9  919  457 (S35E06 L= 85)  1056    251        2 ( 868/ 354- 985/ 556) -69.4 E09-W00 S43-S27
     10 1489  509 (S29W37 L=128)   396     93        2 (1453/ 480-1525/ 546)  59.0 W34-W39 S31-S26
  

    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