Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  862 1478 (N15E21 L=238)  1104    299        1 ( 821/1354- 902/1614) -47.1 E28-E15 N07-N22
      2  312 1450 (N03E60 L=199)  1212    261        1 ( 217/1407- 391/1485)  16.4 E70-E52 N00-N05
      3  395 1415 (N02E50 L=208)   326     42        0 ( 381/1384- 410/1443) -56.3 E52-E48 N00-N03
      4 1476 1338 (N23W19 L=278)   550    122        1 (1452/1268-1490/1394) -76.0 W16-W20 N18-N26
      5 1002 1273 (N06E07 L=251)   515    144        0 ( 929/1226-1060/1323) -16.7 E12-E02 N04-N07
      6  280 1196 (S13E53 L=205)   687    181        1 ( 231/1164- 343/1210)  10.3 E58-E47 S15-S11
      7  193 1104 (S19E60 L=199)  1494    273        4 ( 155/1021- 244/1199) -42.3 E65-E53 S23-S13
      8 1470 1066 (N06W24 L=283)  1030    249        2 (1430/ 933-1507/1189) -53.1 W19-W29 N00-N12
      9  622  873 (S25E21 L=238)   801    284        1 ( 498/ 838- 750/ 923)  32.0 E29-E13 S29-S19
     10 1260  834 (S12W18 L=277)   239    101        0 (1200/ 807-1307/ 863)   0.0 W14-W21 S12-S11
     11  257  792 (S38E50 L=208)  1545    413        3 ( 168/ 742- 376/ 827)  20.0 E60-E38 S41-S33
     12 1633  418 (S22W62 L=320)  1638    199        1 (1565/ 377-1681/ 486)  71.6 W57-W67 S26-S17
  

    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