Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  570 1731 (N41E45 L=295)   452     83        0 ( 527/1711- 604/1748) -10.3 E50-E39 N39-N42
      2  477 1706 (N39E54 L=286)   595     64        0 ( 459/1667- 501/1733) -71.6 E56-E51 N35-N41
      3 1543 1567 (N34W37 L= 16)   370     44        0 (1527/1535-1558/1600)  78.7 W35-W38 N31-N36
      4  324 1547 (N27E60 L=280)   676     46        0 ( 296/1523- 346/1567) -16.7 E65-E55 N25-N28
      5  311  846 (S17E50 L=290)   433    112        0 ( 297/ 792- 335/ 902)  60.3 E50-E46 S20-S13
      6  205  760 (S22E64 L=276)   893     82        0 ( 160/ 738- 252/ 783)   4.1 E71-E57 S23-S20
      7 1714  642 (S22W54 L= 34)   380     67        0 (1693/ 609-1730/ 663) -36.9 W51-W56 S24-S21
      8 1660  528 (S30W54 L= 34)   995     97        3 (1622/ 498-1695/ 559)  76.0 W50-W60 S32-S28
      9  508  513 (S39E39 L=301)   963    129        2 ( 463/ 478- 560/ 554)  20.0 E44-E33 S41-S35
     10 1205  458 (S39W18 L=358)   332     97        0 (1168/ 428-1249/ 502)  56.3 W15-W20 S41-S35
     11 1492  455 (S37W42 L= 22)   421     72        0 (1479/ 415-1509/ 491) -45.0 W39-W45 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