Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  950 1783 (N45E17 L=262)   786    124        2 ( 882/1761-1011/1798)  -4.8 E23-E11 N43-N46
      2 1710 1594 (N42W55 L=334)   572     52        0 (1686/1576-1744/1610)   5.7 W50-W61 N39-N43
      3 1624 1575 (N38W42 L=320)   373     43        0 (1599/1550-1650/1600) -36.9 W39-W44 N36-N40
      4 1165 1538 (N28W04 L=282)   469     87        2 (1123/1511-1200/1583) -51.3 W01-W06 N26-N31
      5  601 1499 (N20E33 L=245)  1112    217        1 ( 523/1428- 671/1571) -28.1 E41-E26 N15-N24
      6 1530 1263 (N15W30 L=308)   508    179        0 (1477/1188-1582/1338) -45.0 W25-W33 N10-N18
      7 1167  825 (S16W11 L=290)   182     93        0 (1134/ 804-1210/ 857) -21.8 W08-W13 S16-S14
      8  150  773 (S26E73 L=206)   590    119        0 ( 131/ 720- 184/ 823)  19.4 E81-E64 S28-S22
      9 1559  776 (S14W38 L=316)   259     69        0 (1539/ 749-1581/ 799) -33.7 W36-W39 S15-S12
     10  819  667 (S29E10 L=269)  1296    392        0 ( 705/ 531- 951/ 796)  49.9 E18-E02 S39-S19
     11 1610  728 (S16W43 L=322)   403     94        2 (1581/ 690-1638/ 767)  63.4 W40-W44 S18-S13
  

    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