Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  466 1665 (N17E57 L=278)   494     97        1 ( 450/1637- 485/1702) -59.0 E60-E55 N15-N20
      2  496 1615 (N15E52 L=283)   521     83        1 ( 459/1587- 530/1648)  45.0 E54-E49 N12-N17
      3 1115 1272 (N09E02 L=334)   476    142        0 (1065/1204-1155/1327) -29.7 E05-W02 N05-N11
      4  477 1204 (S10E37 L=298)   244     40        0 ( 462/1173- 493/1234)  90.0 E38-E36 S11-S08
      5  383 1208 (S12E45 L=290)   379     61        1 ( 355/1183- 414/1228)  26.6 E46-E42 S13-S10
      6 1510 1113 (N11W25 L=  0)   471    122        1 (1490/1045-1527/1180) -60.3 W23-W27 N07-N15
      7 1355 1017 (N01W18 L=353)   677    240        0 (1270/ 904-1394/1113) -30.3 W11-W23 S03-N04
      8  625  833 (S28E19 L=316)   661    252        2 ( 516/ 807- 754/ 857)  13.2 E27-E10 S30-S25
      9 1591  755 (S05W41 L= 16)   243     89        1 (1554/ 723-1613/ 778)  -8.1 W37-W44 S06-S04
     10  856  562 (S38W03 L=338)   318     78        0 ( 818/ 536- 892/ 600)  63.4 W01-W05 S40-S34
     11 1530  398 (S25W55 L= 30)   508     97        0 (1501/ 381-1568/ 423)   8.1 W51-W58 S26-S24
  

    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