Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1412 1629 (N44W33 L=159)  1303    213        1 (1374/1551-1451/1696)  50.2 W28-W40 N37-N49
      2 1020 1628 (N44E01 L=125)  1078    187        3 ( 972/1544-1058/1694) -85.2 E06-W01 N37-N49
      3  356 1544 (N34E57 L= 69)   569     90        1 ( 335/1514- 389/1566)  26.6 E59-E53 N32-N35
      4  572 1460 (N30E34 L= 92)   396    109        0 ( 556/1409- 592/1509) -74.1 E35-E32 N26-N33
      5 1605 1449 (N30W44 L=170)   296     49        1 (1575/1430-1630/1469) -14.0 W41-W45 N29-N31
      6 1446 1436 (N30W30 L=156)   325     79        2 (1419/1398-1467/1465) -53.1 W28-W31 N27-N31
      7 1134 1330 (N22W07 L=133)   613    199        2 (1049/1262-1227/1397) -36.0 W01-W12 N18-N26
      8  650 1194 (N13E24 L=102)   382    114        0 ( 619/1144- 685/1252) -50.2 E27-E21 N09-N16
      9 1123 1214 (N15W06 L=132)   255    111        2 (1074/1196-1190/1229)   7.1 W02-W10 N14-N16
     10 1701 1193 (N14W47 L=173)   350     70        1 (1673/1164-1724/1213)   0.0 W44-W49 N11-N14
     11 1096 1119 (N09W04 L=130)   581    242        0 (1008/1074-1225/1189) -28.3 E01-W12 N06-N13
     12 1468  881 (S05W29 L=155)   352    111        1 (1449/ 827-1489/ 934) -80.5 W27-W30 S08-S02
     13  383  855 (S09E44 L= 82)   333    108        0 ( 343/ 828- 430/ 885) -23.2 E47-E39 S10-S07
     14 1199  722 (S15W12 L=138)   229     93        0 (1181/ 673-1211/ 771)  90.0 W10-W12 S18-S12
     15  339  676 (S21E50 L= 76)   474    113        0 ( 306/ 654- 382/ 713) -26.6 E52-E45 S22-S18
     16  564  388 (S41E38 L= 88)  2326    423        3 ( 419/ 338- 661/ 442)   9.5 E52-E27 S45-S36
     17 1455  378 (S40W38 L=164)  2127    347        2 (1307/ 337-1598/ 437)   9.5 W25-W49 S43-S35
  

    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