Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  564 1731 (N41E46 L=294)   489     56        0 ( 528/1714- 597/1746)  -6.3 E50-E41 N39-N42
      2  472 1704 (N39E55 L=285)   594     57        0 ( 453/1668- 495/1733) -78.7 E57-E52 N36-N41
      3 1536 1567 (N34W36 L= 16)   353     57        0 (1519/1535-1550/1605)  63.4 W34-W38 N31-N37
      4  319 1548 (N27E61 L=279)   519     34        1 ( 297/1527- 339/1565) -20.6 E65-E57 N25-N28
      5  309  850 (S17E50 L=290)   520    105        0 ( 293/ 797- 334/ 907)  56.3 E51-E46 S20-S13
      6  203  762 (S22E64 L=276)   852     79        0 ( 161/ 740- 248/ 782)   0.0 E71-E58 S23-S20
      7 1394  650 (S25W27 L=  8)   217     65        0 (1365/ 632-1430/ 662)  36.9 W25-W29 S26-S23
      8 1711  644 (S22W54 L= 34)   336     30        0 (1691/ 620-1728/ 662) -33.7 W51-W56 S23-S21
      9  508  516 (S38E39 L=301)   906    132        1 ( 466/ 483- 561/ 558)  20.0 E44-E33 S40-S35
     10 1656  530 (S30W54 L= 34)   723     84        3 (1621/ 502-1693/ 558)  63.4 W50-W59 S32-S28
     11 1203  460 (S39W17 L=358)   303     89        0 (1165/ 432-1244/ 500)  38.7 W14-W19 S41-S35
     12 1487  456 (S37W41 L= 22)   310     64        0 (1476/ 424-1498/ 496) -51.3 W39-W44 S39-S34
  

    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