Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1255 1437 (N29W18 L=299)   488    150        1 (1218/1369-1285/1499)  83.7 W15-W20 N24-N33
      2 1427 1381 (N25W30 L=311)   242     82        0 (1389/1352-1453/1405) -45.0 W27-W31 N22-N26
      3 1341 1350 (N23W23 L=304)   400     60        0 (1305/1331-1385/1369) -18.4 W20-W26 N21-N24
      4  811 1340 (N25E13 L=268)   220     86        1 ( 782/1313- 834/1364)  36.9 E15-E11 N23-N26
      5  420 1221 (N18E42 L=239)   265     73        0 ( 400/1187- 437/1259)  90.0 E43-E40 N15-N20
      6  930 1151 (N13E05 L=275)   266     83        1 ( 895/1119- 957/1183)  53.1 E07-E03 N10-N14
      7 1243  824 (S09W13 L=294)   513    118        3 (1220/ 777-1270/ 883)  63.4 W11-W14 S11-S05
      8 1748  855 (S11W51 L=332)   548     39        1 (1724/ 827-1776/ 882)  90.0 W48-W54 S12-S08
      9 1585  796 (S13W37 L=318)   289    119        0 (1565/ 735-1601/ 844)  81.9 W35-W38 S16-S09
     10 1640  748 (S17W42 L=323)   730    116        4 (1603/ 707-1693/ 792)  23.2 W38-W46 S18-S13
     11 1750  661 (S24W55 L=336)  2271    320        3 (1692/ 588-1800/ 764) -49.8 W48-W64 S29-S16
     12  280  598 (S22E64 L=217)  1064    181        1 ( 258/ 524- 328/ 670)  79.7 E68-E59 S27-S16
  

    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