Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1513 1714 (N31W60 L=234)   856    134        1 (1456/1691-1574/1734) -14.0 W53-W65 N27-N32
      2 1109 1656 (N39W25 L=199)   923    233        1 (1000/1614-1209/1694)  -6.3 W15-W33 N37-N39
      3  351 1282 (N33E42 L=132)   933    204        1 ( 263/1264- 450/1303) -16.7 E52-E32 N30-N36
      4  522 1159 (N22E28 L=146)   248     76        0 ( 484/1140- 561/1172)   0.0 E31-E25 N20-N22
      5 1654 1012 (S12W38 L=213)   366    130        2 (1628/ 948-1680/1068)  45.0 W35-W41 S15-S09
      6  176  970 (N18E63 L=111)  1703    242        1 (  86/ 918- 253/1043)   5.9 E82-E52 N14-N20
      7 1772 1004 (S16W48 L=223)   433     63        1 (1758/ 969-1785/1040)  36.9 W46-W50 S17-S13
      8 1834  902 (S24W54 L=228)  1008    135        2 (1795/ 857-1861/ 969)  45.0 W49-W57 S26-S19
      9 1087  939 (S03W02 L=176)   216     76        0 (1042/ 927-1129/ 953) -33.7 W00-W04 S03-S01
     10 1387  727 (S22W15 L=189)   250     74        0 (1374/ 682-1401/ 763)  53.1 W13-W16 S24-S20
     11 1191  620 (S24W00 L=174)   711    253        1 (1136/ 522-1248/ 732)  42.3 E05-W06 S29-S18
     12  994  596 (S21E12 L=162)   802    246        1 ( 957/ 485-1039/ 712)  47.5 E17-E06 S27-S14
     13  708  612 (S13E29 L=145)   265     68        0 ( 674/ 580- 739/ 643)  31.0 E31-E26 S14-S11
     14 1253  532 (S31W02 L=176)   351    128        0 (1195/ 505-1299/ 570)   7.1 E02-W06 S32-S29
     15  671  398 (S26E41 L=133)   883    223        0 ( 634/ 303- 697/ 512)  40.9 E49-E34 S31-S18
     16 1529  404 (S50W23 L=197)   789    119        0 (1466/ 383-1584/ 436)  -8.7 W16-W29 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