Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  746 1697 (N40E35 L= 73)   618    109        0 ( 715/1665- 785/1735) -18.4 E39-E30 N37-N42
      2  451 1638 (N31E56 L= 52)   614    186        0 ( 436/1597- 470/1691) -45.0 E60-E51 N27-N34
      3 1291 1631 (N44W11 L=119)   435     59        0 (1275/1598-1307/1666)  80.5 W09-W12 N41-N47
      4  766 1590 (N32E28 L= 80)   766    186        3 ( 741/1521- 800/1658) -78.7 E29-E26 N27-N37
      5 1630 1573 (N46W44 L=152)   908    202        0 (1587/1519-1660/1620)  39.8 W37-W49 N40-N50
      6  619 1478 (N22E35 L= 73)   902    271        1 ( 503/1450- 721/1503)  18.4 E43-E27 N20-N25
      7  750 1365 (N17E22 L= 86)   444    212        0 ( 682/1319- 784/1436)  82.9 E26-E19 N13-N21
      8  349 1249 (N05E49 L= 59)   431     87        1 ( 311/1227- 388/1274)  26.6 E52-E46 N03-N06
      9  901 1187 (N08E10 L= 98)   189     32        1 ( 883/1162- 915/1211)  76.0 E10-E08 N06-N10
     10 1390 1167 (N14W21 L=129)   444    162        1 (1337/1127-1421/1211)  50.2 W17-W23 N10-N16
     11 1329 1107 (N09W18 L=126)   260     52        0 (1304/1087-1359/1126)   0.0 W15-W19 N08-N10
     12  383  971 (S12E42 L= 66)   306     94        0 ( 371/ 924- 397/1023) -80.5 E43-E40 S14-S08
     13 1800  891 (N03W57 L=165)   426     83        0 (1787/ 849-1814/ 935) -78.7 W55-W58 N00-N05
     14  296  801 (S24E49 L= 59)   866    411        1 ( 194/ 731- 386/ 883)  31.0 E61-E41 S29-S17
     15  990  838 (S11W00 L=108)   227     77        0 ( 949/ 819-1033/ 861) -11.3 E02-W03 S12-S10
     16 1392  711 (S14W28 L=136)   641    161        5 (1350/ 653-1416/ 757) -39.8 W24-W30 S17-S11
     17 1749  686 (S10W58 L=166)   483     90        0 (1726/ 637-1779/ 727)  80.5 W55-W60 S13-S07
  

    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