Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1034 1748 (N52W29 L=290)   580     86        1 ( 997/1715-1065/1780)   0.0 W25-W34 N48-N53
      2 1462 1717 (N34W59 L=320)   414     60        0 (1436/1705-1497/1734) -45.0 W56-W62 N32-N35
      3  923 1680 (N50W15 L=275)   346     53        0 ( 892/1658- 950/1697) -36.9 W12-W17 N48-N52
      4 1671 1623 (N22W71 L=331)   569     47        2 (1646/1603-1693/1646)  33.7 W67-W75 N20-N22
      5  542 1356 (N38E26 L=235)  1177    230        4 ( 451/1282- 630/1396) -29.4 E34-E18 N31-N41
      6 1048 1332 (N24W09 L=270)  1137    382        3 ( 908/1303-1228/1359) -19.3 W00-W20 N20-N27
      7 1760 1142 (S06W51 L=311)  1121    146        3 (1726/1075-1792/1204)  40.6 W46-W55 S09-S02
      8  567 1072 (N20E27 L=233)   983    239        1 ( 475/ 987- 656/1135) -50.2 E33-E22 N13-N25
      9  177  987 (N21E66 L=194)  2392    272        4 (  97/ 902- 249/1096) -22.5 E86-E57 N15-N27
     10 1665  621 (S35W35 L=295)  6700    821       10 (1463/ 382-1856/ 973)  33.1 W14-W57 S49-S18
     11 1280  900 (S06W12 L=272)   287     59        0 (1260/ 865-1305/ 941) -90.0 W10-W12 S08-S03
     12 1221  860 (S07W07 L=267)   270    125        0 (1166/ 847-1281/ 871) -33.7 W03-W10 S09-S05
     13 1251  657 (S20W04 L=265)   650    251        0 (1190/ 569-1282/ 763)  51.3 E01-W07 S24-S13
     14 1344  311 (S49W02 L=262)   413     50        2 (1322/ 286-1365/ 343)  68.2 W00-W03 S51-S45
     15 1258  262 (S51E07 L=253)   686     73        2 (1222/ 237-1297/ 285) -39.8 E09-E03 S53-S47
     16  932  238 (S41E33 L=227)   556     73        1 ( 894/ 215- 973/ 265)   6.3 E37-E28 S42-S39
     17 1152  211 (S52E19 L=241)   572     74        1 (1120/ 185-1192/ 233)   0.0 E23-E14 S54-S50
  

    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