Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  956 1567 (N38W13 L=162)   625    181        1 ( 909/1487-1013/1630)  28.6 W08-W19 N34-N41
      2 1244 1548 (N28W30 L=179)   250     47        0 (1227/1517-1261/1576)  76.0 W28-W31 N26-N30
      3 1371 1528 (N23W37 L=186)  1131    197        4 (1288/1495-1457/1557) -31.0 W31-W41 N19-N26
      4  416 1184 (N30E35 L=114)  1055    373        3 ( 280/1144- 548/1236) -22.2 E46-E24 N25-N36
      5  719 1095 (N17E16 L=133)   677    150        2 ( 658/1074- 797/1123) -32.0 E19-E10 N14-N20
      6 1373 1073 (S01W21 L=170)  1018    231        2 (1273/1040-1491/1108) -24.8 W14-W27 S06-N01
      7 1333  995 (S04W16 L=165)   347    117        1 (1279/ 955-1392/1013)  -7.1 W12-W20 S05-S03
      8 1543  907 (S15W27 L=176)   540     96        1 (1499/ 863-1587/ 935)  -8.1 W23-W30 S17-S13
      9 1261  760 (S15W06 L=155)   734    186        1 (1231/ 679-1297/ 859)  63.4 W04-W09 S20-S10
     10  935  611 (S15E16 L=133)   681    152        1 ( 840/ 576-1000/ 637) -48.8 E21-E13 S19-S11
     11  630  387 (S21E45 L=105)  1047    224        4 ( 563/ 310- 697/ 465) -85.2 E46-E42 S27-S15
     12  917  181 (S45E41 L=108)  1154    176        2 ( 866/ 135- 951/ 214) -18.4 E47-E36 S49-S41
  

    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