Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1116 1481 (N32W05 L=332)   436    121        0 (1063/1439-1153/1524) -45.0 W01-W08 N29-N35
      2 1405 1449 (N31W27 L=354)   535    208        0 (1322/1389-1481/1485) -22.6 W20-W32 N26-N33
      3 1699 1413 (N28W53 L= 20)   865    191        3 (1664/1348-1740/1484)  71.6 W49-W58 N23-N32
      4 1375 1398 (N27W24 L=351)   263    118        0 (1322/1372-1423/1419) -23.2 W20-W27 N25-N28
      5  417 1193 (N11E42 L=285)   249     98        2 ( 394/1167- 453/1229) -38.7 E43-E38 N09-N13
      6 1712 1181 (N13W48 L= 15)   551     82        3 (1672/1164-1759/1202)   0.0 W44-W53 N12-N14
      7 1618  831 (S08W40 L=  7)   480    174        2 (1588/ 774-1667/ 903) -48.8 W37-W45 S11-S03
      8 1349  841 (S08W21 L=348)   271     79        2 (1317/ 822-1383/ 874) -11.3 W18-W23 S08-S05
      9 1741  776 (S12W52 L= 19)   435     77        1 (1723/ 733-1756/ 810) -68.2 W50-W54 S14-S09
     10 1526  710 (S16W35 L=  2)   422    114        2 (1496/ 674-1560/ 746) -38.7 W32-W37 S18-S13
     11 1211  683 (S18W13 L=340)   258     83        0 (1179/ 655-1243/ 714)  45.0 W10-W14 S19-S15
     12 1398  582 (S25W27 L=354)  2198    522       10 (1179/ 445-1537/ 695)  37.3 W12-W36 S34-S16
     13 1339  675 (S18W22 L=349)   268    110        0 (1287/ 658-1383/ 692)   0.0 W17-W24 S19-S17
     14  344  619 (S25E52 L=275)   841    165        1 ( 303/ 558- 381/ 689)  35.5 E60-E46 S30-S20
     15  523  528 (S31E37 L=290)  1667    273        3 ( 400/ 488- 651/ 584)  15.1 E51-E25 S34-S26
     16 1416  455 (S34W32 L=359)   593    195        0 (1324/ 435-1504/ 481)  12.1 W24-W39 S35-S31
  

    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