Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  733 1600 (N28E26 L=129)   845    276        0 ( 676/1555- 843/1674)   4.8 E30-E17 N25-N33
      2 1740 1379 (N25W48 L=202)   835     89        2 (1701/1336-1778/1426) -51.3 W44-W51 N22-N27
      3  317 1366 (N11E52 L=103)   446     84        1 ( 289/1325- 351/1404)  68.2 E53-E49 N08-N13
      4 1234 1352 (N16W10 L=164)   242     49        1 (1203/1330-1256/1375)   0.0 W07-W11 N15-N17
      5 1330 1318 (N15W16 L=170)   297     59        1 (1303/1285-1356/1341)  71.6 W14-W17 N13-N16
      6 1670 1312 (N19W40 L=195)   269     53        0 (1652/1288-1691/1341) -56.3 W39-W42 N17-N20
      7 1126 1209 (N06W04 L=159)  1221    394        2 ( 979/1058-1199/1337) -45.0 E05-W09 S01-N12
      8 1676 1215 (N13W41 L=195)   265     91        0 (1663/1176-1695/1248)  59.0 W39-W42 N10-N15
      9  563 1131 (S02E29 L=125)   525    159        1 ( 503/1086- 630/1165) -12.5 E33-E24 S04-N00
     10 1027 1044 (S04W00 L=154)   220     72        0 ( 992/1026-1074/1060) -14.0 E02-W02 S04-S03
     11  870  973 (S09E09 L=146)   272    118        0 ( 827/ 935- 924/1004)  38.7 E11-E05 S12-S07
     12  198  931 (S15E59 L= 95)   404     90        0 ( 180/ 888- 221/ 967)  38.7 E61-E56 S17-S13
     13  277  892 (S18E51 L=103)   652    143        1 ( 250/ 825- 298/ 956) -82.9 E53-E49 S21-S13
     14  139  873 (S19E68 L= 86)  1277    168        4 ( 108/ 819- 169/ 947) -74.1 E74-E64 S21-S14
     15 1419  894 (S09W25 L=180)   214     13        0 (1398/ 876-1439/ 912) -18.4 W23-W26 S09-S07
     16 1113  785 (S18W08 L=162)  1480    709        1 ( 855/ 719-1415/ 861)   8.4 E08-W26 S23-S13
     17  380  652 (S32E45 L=109)   864    170        3 ( 308/ 619- 460/ 698)  14.0 E53-E37 S34-S28
     18 1413  462 (S35W36 L=190)   803    133        1 (1349/ 441-1490/ 498)   0.0 W30-W43 S36-S33
     19 1415  360 (S42W42 L=197)   747    131        1 (1378/ 327-1480/ 408)  -4.4 W37-W50 S44-S38
  

    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