Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  847 1729 (N35E34 L=348)   983    193        2 ( 816/1670- 876/1795) -59.0 E38-E30 N30-N40
      2 1387 1394 (N27W13 L= 35)  1710    426        3 (1207/1270-1564/1535) -13.6 E02-W27 N23-N32
      3  655 1383 (N08E31 L=351)   923    267        1 ( 602/1317- 756/1481) -19.7 E37-E23 N05-N13
      4 1730 1382 (N38W42 L= 63)   621    126        1 (1691/1355-1775/1403)   0.0 W37-W46 N35-N39
      5 1031 1291 (N12E06 L= 16)   488    104        0 ( 992/1241-1065/1342) -39.8 E09-E03 N09-N14
      6 1148 1145 (N06W04 L= 26)   185     83        0 (1139/1102-1158/1189) -53.1 W02-W05 N04-N08
      7 1110 1080 (N02W03 L= 25)   200     80        0 (1064/1069-1141/1100)  21.8 W00-W05 N00-N03
      8  947  803 (S18W01 L= 23)   193     59        0 ( 916/ 778- 969/ 837) -33.7 E01-W02 S18-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