Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  582 1801 (N37E67 L=148)   682     62        0 ( 561/1774- 611/1828)  -9.0 E79-E59 N35-N38
      2 1273 1603 (N40W02 L=218)  1070    206        3 (1212/1515-1317/1676) -48.0 E04-W06 N33-N44
      3 1207 1398 (N24W02 L=218)  1145    614        2 (1062/1143-1309/1592) -42.6 E12-W13 N10-N34
      4  760 1546 (N23E30 L=186)   278     86        0 ( 729/1512- 785/1587) -23.2 E34-E27 N21-N24
      5  697 1444 (N15E30 L=186)   307    100        1 ( 673/1405- 720/1494) -33.7 E33-E27 N12-N17
      6  824 1395 (N15E21 L=195)   343    117        1 ( 779/1362- 869/1430)  -9.5 E23-E17 N13-N16
      7 1104 1359 (N19E03 L=213)   230    117        0 (1085/1305-1119/1414) -81.9 E04-E01 N15-N22
      8 1561 1277 (N25W29 L=245)   792    267        1 (1521/1205-1652/1327)   0.0 W25-W36 N20-N28
      9  143 1203 (S10E75 L=141)  1339    112        2 ( 111/1146- 172/1247)  35.0 E81-E70 S13-S06
     10  232 1089 (S15E58 L=158)  1963    347        5 ( 197/ 959- 274/1199)  90.0 E63-E53 S23-S08
     11 1795 1065 (N18W53 L=269)   937    167        2 (1780/ 975-1813/1138) -78.7 W50-W54 N12-N22
     12  765  962 (S12E14 L=202)   220     93        1 ( 733/ 940- 797/ 987) -11.3 E16-E11 S12-S09
     13 1221  889 (S06W14 L=230)   308    121        0 (1184/ 840-1246/ 936) -33.7 W11-W17 S08-S03
     14  971  760 (S19W03 L=219)   292     87        0 ( 924/ 748-1020/ 775)  31.0 W00-W05 S20-S17
     15  328  698 (S38E47 L=169)  3008    358        6 ( 170/ 656- 473/ 747)   7.8 E83-E31 S43-S32
     16  703  394 (S51E07 L=208)   525    118        0 ( 660/ 374- 759/ 424)   0.0 E12-E01 S52-S49
     17 1507  330 (S32W62 L=278)   710     46        0 (1483/ 308-1536/ 357)  -5.7 W56-W66 S33-S30
  

    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