Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  947 1731 (N37E29 L=300)  1128    224        3 ( 902/1635- 999/1818) -30.1 E39-E20 N31-N43
      2  646 1695 (N26E48 L=281)   788    166        2 ( 610/1620- 691/1761) -26.6 E56-E40 N22-N30
      3 1504 1701 (N56W18 L=346)  1672    221        3 (1425/1657-1577/1732)  30.1 W09-W28 N48-N61
      4  176 1000 (S25E64 L=265)  1450    208        2 ( 123/ 915- 238/1064)  40.2 E74-E55 S31-S20
      5  518  753 (S33E26 L=302)   685    191        1 ( 433/ 725- 623/ 773)  25.0 E33-E18 S35-S28
      6 1037  606 (S27W13 L=341)   553    124        0 (1008/ 545-1074/ 659)  90.0 W11-W14 S31-S23
      7 1259  457 (S30W33 L=  2)   271    134        0 (1220/ 430-1308/ 510) -12.5 W28-W37 S30-S27
      8  355  466 (S59E53 L=276)   500     33        1 ( 332/ 448- 378/ 488)  22.6 E59-E46 S61-S56
      9  420  429 (S61E41 L=287)   980    118        1 ( 390/ 396- 467/ 481) -56.3 E48-E33 S63-S55
     10  516  375 (S62E27 L=302)   796     62        1 ( 489/ 341- 539/ 412) -81.9 E31-E22 S65-S58
     11  683  332 (S59E03 L=326)   340     15        0 ( 667/ 311- 697/ 356) -36.9 E05-W00 S61-S56
     12  880  223 (S61W26 L=354)   466     54        0 ( 861/ 195- 900/ 255) -20.0 W20-W31 S63-S58
  

    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