Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  605 1564 (N37E35 L=240)   609    184        1 ( 555/1525- 670/1607) -26.6 E40-E28 N34-N40
      2 1051 1489 (N33W01 L=275)   274     59        2 (1018/1470-1086/1507)  36.9 E01-W03 N31-N34
      3  291 1291 (N17E55 L=219)   522     89        1 ( 267/1266- 325/1313)  36.9 E58-E51 N15-N18
      4  756 1220 (N15E18 L=257)   538    232        1 ( 725/1139- 815/1299) -78.7 E19-E13 N09-N19
      5  899 1185 (N13E08 L=266)   421    172        1 ( 826/1140- 966/1258)  41.2 E12-E04 N10-N17
      6  155 1160 (N08E70 L=204)  1622    259        4 ( 128/1092- 188/1231) -82.9 E76-E65 N03-N11
      7  816  859 (S07E13 L=262)   346    156        0 ( 770/ 814- 874/ 881)   0.0 E15-E09 S10-S05
      8 1205  751 (S14W12 L=286)   427    118        1 (1176/ 697-1234/ 813)  69.4 W10-W13 S17-S09
      9  883  783 (S12E08 L=266)   209     49        0 ( 853/ 759- 908/ 802) -14.0 E10-E06 S13-S10
     10  237  676 (S22E61 L=213)   498    117        1 ( 208/ 656- 260/ 694) -26.6 E64-E57 S23-S20
     11  554  584 (S26E33 L=242)   345    123        2 ( 539/ 536- 581/ 628)  63.4 E34-E30 S29-S23
     12 1094  564 (S26W06 L=280)   461    135        0 (1023/ 541-1151/ 604)  18.4 W00-W09 S27-S23
  

    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