Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1149 1583 (N40W06 L=340)   488     68        0 (1103/1561-1188/1601)   0.0 W01-W08 N38-N41
      2  672 1476 (N29E29 L=305)   734    213        0 ( 592/1430- 787/1512) -13.2 E36-E19 N26-N30
      3 1446 1426 (N30W28 L=  2)   374     89        0 (1411/1401-1485/1450)  -9.5 W25-W31 N28-N31
      4  501 1307 (N16E38 L=296)  1220    268        2 ( 427/1230- 595/1408) -27.9 E46-E29 N11-N22
      5 1706 1347 (N26W51 L= 25)   361     83        0 (1686/1305-1722/1383) -78.7 W49-W52 N23-N28
      6  977 1303 (N19E05 L=329)   602    206        0 ( 885/1247-1049/1349) -24.4 E11-W00 N16-N21
      7  796  723 (S18E13 L=321)   303     67        1 ( 761/ 702- 836/ 743) -11.3 E15-E10 S19-S16
      8 1457  696 (S16W31 L=  5)   355     68        1 (1431/ 672-1482/ 727)  53.1 W28-W32 S17-S13
      9 1707  652 (S18W54 L= 28)   409     51        0 (1688/ 628-1731/ 676)  56.3 W51-W55 S19-S16
     10  298  542 (S35E63 L=271)   786     89        0 ( 252/ 522- 343/ 569)   9.9 E77-E54 S37-S32
  

    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