Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1529 1531 (N29W47 L= 62)   840    234        2 (1461/1458-1605/1570)   0.0 W38-W52 N25-N31
      2 1153 1445 (N30W15 L= 30)   627    258        1 (1126/1344-1186/1538)  85.6 W12-W18 N23-N36
      3  425 1294 (N29E41 L=334)  1894    379        5 ( 391/1140- 498/1460) -59.4 E46-E33 N18-N40
      4 1648 1372 (N17W50 L= 65)   568    135        3 (1609/1305-1687/1420) -82.9 W47-W51 N12-N20
      5 1530  986 (S04W31 L= 46)   435    149        1 (1501/ 918-1568/1044)  54.5 W29-W35 S07-N00
      6 1303  967 (S01W16 L= 31)   217     64        1 (1269/ 954-1338/ 981) -14.0 W14-W18 S02-N00
      7 1747  859 (S16W48 L= 63)   568     68        1 (1711/ 834-1772/ 884)   8.1 W44-W51 S17-S14
      8  723  577 (S18E26 L=349)   261     66        0 ( 692/ 556- 757/ 600)  11.3 E28-E23 S19-S16
      9  624  559 (S18E34 L=341)   356     57        1 ( 601/ 537- 655/ 587) -45.0 E35-E31 S20-S16
     10  568  508 (S21E40 L=335)   284     97        0 ( 544/ 479- 610/ 525)   0.0 E41-E36 S23-S20
     11  952  486 (S27E13 L=  2)   307    113        1 ( 902/ 468-1013/ 497)  -7.1 E16-E08 S28-S26
  

    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