Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1083 1906 (N58W58 L=315)   721     73        1 (1049/1880-1117/1927) -76.0 W51-W65 N55-N60
      2 1003 1458 (N30W12 L=268)   210     64        0 ( 969/1443-1035/1471) -14.0 W09-W13 N29-N31
      3 1619  896 (S18W32 L=289)   293     74        0 (1601/ 856-1634/ 938)  78.7 W31-W33 S20-S15
      4 1703  841 (S24W38 L=294)   697    183        0 (1638/ 786-1766/ 901)  13.0 W32-W45 S25-S22
      5 1492  664 (S28W19 L=275)   296     94        0 (1477/ 619-1505/ 713)  71.6 W17-W20 S31-S25
      6 1428  628 (S29W14 L=270)   487    127        0 (1388/ 575-1458/ 681)  45.0 W09-W16 S31-S25
      7 1559  623 (S33W23 L=279)   451     86        0 (1511/ 608-1603/ 642) -20.6 W19-W27 S35-S31
      8 1234  568 (S26E00 L=256)   486    121        0 (1208/ 513-1260/ 618)  50.2 E03-W02 S29-S23
      9 1769  525 (S51W50 L=307)   573     57        0 (1735/ 512-1802/ 538) -11.3 W43-W59 S53-S48
     10  589  375 (S22E49 L=207)   507     70        1 ( 561/ 344- 618/ 414)  51.3 E51-E46 S23-S18
  

    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