Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  472 1671 (N17E59 L=  6)   451     88        1 ( 447/1645- 496/1694) -23.2 E61-E54 N15-N18
      2  418 1602 (N12E59 L=  6)   447     63        1 ( 398/1574- 440/1628) -45.0 E60-E56 N10-N13
      3  606 1563 (N13E42 L= 23)   261     29        0 ( 590/1539- 620/1584)  90.0 E43-E40 N11-N14
      4 1434 1291 (N20W16 L= 81)   911    273        0 (1411/1147-1460/1413) -65.0 W12-W20 N12-N28
      5  519 1303 (S04E38 L= 27)   297     87        0 ( 505/1263- 534/1347) -59.0 E39-E35 S06-S01
      6 1736 1278 (N30W39 L=104)   481     52        1 (1708/1248-1757/1306) -33.7 W36-W40 N28-N32
      7 1611 1129 (N17W32 L= 96)   798    151        1 (1583/1045-1630/1207) -56.3 W27-W34 N11-N20
      8  261  927 (S31E49 L= 15)  1639    430        3 ( 140/ 905- 378/ 956)   7.1 E62-E38 S33-S28
      9 1447  831 (S05W29 L= 94)   215     35        0 (1426/ 805-1470/ 853)  76.0 W28-W30 S06-S02
     10  877  798 (S23E02 L= 63)   407    102        0 ( 836/ 756- 924/ 830)  71.6 E04-W00 S26-S19
     11  522  683 (S40E24 L= 41)   757    166        1 ( 453/ 659- 599/ 705)  26.6 E29-E17 S42-S36
     12 1082  517 (S33W20 L= 85)  1116    195        0 ( 993/ 457-1177/ 592)  69.0 W16-W23 S38-S25
     13 1424  512 (S22W41 L=106)   747    139        1 (1350/ 472-1483/ 561)  63.4 W37-W43 S25-S17
     14 1651  519 (S14W58 L=123)   439    158        2 (1605/ 492-1692/ 548)  11.3 W52-W62 S15-S12
     15 1545  440 (S22W54 L=119)   359     57        0 (1514/ 418-1576/ 468)  63.4 W52-W55 S23-S19
     16 1387  372 (S31W47 L=112)   606    107        1 (1325/ 354-1432/ 389)  39.8 W43-W50 S33-S28
  

    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