Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  941 1735 (N37E30 L=300)  1141    221        3 ( 903/1640- 990/1826) -33.7 E39-E21 N31-N43
      2  644 1687 (N26E47 L=282)   871    168        4 ( 608/1620- 687/1768) -25.2 E57-E40 N21-N30
      3 1503 1704 (N57W18 L=347)  1745    177        3 (1426/1659-1588/1734)  17.7 W08-W31 N50-N61
      4  177 1005 (S25E64 L=266)  1338    208        2 ( 131/ 916- 249/1065)  40.2 E73-E54 S31-S20
      5  511  756 (S33E27 L=303)   682    190        2 ( 428/ 727- 614/ 778)  26.6 E34-E19 S35-S28
      6 1028  610 (S27W12 L=341)   563    119        0 ( 998/ 549-1064/ 666) -82.9 W10-W13 S31-S23
      7 1301  463 (S28W36 L=  5)   509    259        0 (1212/ 432-1397/ 517)  12.1 W27-W41 S31-S24
      8  414  431 (S61E43 L=287)   999    125        1 ( 377/ 400- 459/ 479) -51.3 E50-E34 S63-S56
      9  511  378 (S62E27 L=302)   792     70        0 ( 487/ 342- 533/ 413) -90.0 E31-E22 S65-S58
     10  679  334 (S59E04 L=326)   343     37        0 ( 664/ 309- 696/ 359) -53.1 E05-E01 S61-S56
     11  877  228 (S61W25 L=354)   507     70        0 ( 855/ 197- 906/ 263) -17.1 W19-W32 S62-S57
  

    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