Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1560 1702 (N32W60 L= 25)   801    153        1 (1521/1680-1605/1722) -33.7 W56-W64 N29-N33
      2 1026 1590 (N41W16 L=340)   407     75        0 (1007/1547-1045/1626)  59.0 W13-W18 N38-N43
      3 1148 1576 (N37W23 L=348)   256     63        0 (1109/1562-1182/1589) -18.4 W20-W26 N36-N38
      4 1021 1425 (N31W09 L=334)   342     92        0 ( 969/1408-1074/1444) -33.7 W05-W12 N28-N32
      5 1667 1386 (N12W52 L= 16)   360     77        3 (1645/1355-1686/1416)  76.0 W49-W53 N09-N13
      6 1309 1246 (N13W22 L=347)   249     47        1 (1281/1217-1335/1268) -76.0 W20-W23 N11-N15
      7 1011 1186 (N17W03 L=327)   263    145        0 ( 940/1167-1079/1200) -23.2 E01-W06 N14-N18
      8 1335 1050 (N02W19 L=343)   363     67        1 (1305/1006-1365/1090) -78.7 W17-W19 S01-N04
      9  280  849 (N09E55 L=270)   390     89        0 ( 257/ 820- 314/ 894) -59.0 E56-E51 N07-N12
     10 1445  725 (S20W19 L=344)   417    114        0 (1407/ 673-1489/ 781)  35.5 W16-W23 S22-S17
     11 1215  410 (S36E02 L=323)   256     79        0 (1201/ 384-1243/ 445)  80.5 E03-W00 S38-S32
  

    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