Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1502 1715 (N31W59 L=233)  1008    207        1 (1391/1692-1577/1737) -20.6 W49-W66 N28-N34
      2 1115 1659 (N39W25 L=199)   955    225        0 (1001/1615-1208/1699)  -3.0 W15-W34 N37-N39
      3 1312 1445 (N20W28 L=202)   261     85        1 (1272/1421-1348/1468) -45.0 W26-W30 N18-N22
      4  349 1285 (N33E42 L=132)   947    178        3 ( 265/1263- 440/1310) -14.7 E52-E33 N30-N36
      5  527 1161 (N22E28 L=146)   273     86        0 ( 486/1143- 571/1174)   0.0 E30-E24 N20-N22
      6 1666 1034 (S11W40 L=214)   480    165        2 (1635/ 955-1698/1106)  48.8 W36-W43 S15-S07
      7  176  970 (N18E63 L=111)  2077    247        1 (  86/ 917- 260/1049)   5.7 E82-E51 N14-N21
      8 1775 1004 (S16W49 L=223)   433     56        0 (1763/ 969-1789/1041)  56.3 W47-W50 S17-S13
      9 1836  911 (S23W55 L=228)  1297    171        3 (1795/ 856-1865/1005)  63.4 W49-W58 S26-S17
     10 1095  943 (S03W02 L=176)   260     97        0 (1044/ 931-1142/ 955) -26.6 W00-W04 S03-S01
     11 1393  729 (S22W15 L=189)   271     60        0 (1378/ 690-1408/ 764)  63.4 W13-W16 S24-S20
     12 1197  624 (S24W01 L=175)   765    243        1 (1143/ 520-1254/ 736)  47.7 E04-W06 S29-S18
     13 1006  613 (S20E11 L=163)   949    266        2 ( 965/ 482-1066/ 719)  45.0 E16-E04 S27-S14
     14  685  469 (S21E36 L=138)  1360    383        2 ( 638/ 302- 745/ 645)  39.8 E49-E25 S32-S11
     15 1257  534 (S31W03 L=176)   379    114        0 (1202/ 505-1306/ 572)  15.9 E01-W06 S32-S29
     16 1563  421 (S50W27 L=201)  1256    205        1 (1465/ 384-1658/ 477)   0.0 W16-W38 S52-S47
  

    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