Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  528 1649 (N43E50 L= 42)   461     44        0 ( 506/1625- 550/1678) -16.7 E55-E45 N41-N45
      2 1410 1571 (N40W31 L=122)   388     64        0 (1378/1547-1439/1593) -31.0 W28-W33 N37-N41
      3  297 1510 (N31E69 L= 22)  1047    158        0 ( 277/1441- 332/1574) -23.2 E79-E58 N26-N35
      4 1679 1448 (N31W53 L=144)   729    148        1 (1636/1378-1717/1489) -63.4 W49-W56 N25-N33
      5 1214 1369 (N25W12 L=104)   361    114        1 (1175/1329-1266/1401)  29.7 W09-W16 N22-N27
      6 1490 1236 (N17W31 L=122)   871    337        1 (1383/1136-1612/1327) -38.7 W23-W39 N10-N22
      7 1855 1262 (N18W68 L=160)  1272     98        1 (1829/1209-1875/1310)  71.6 W65-W72 N14-N20
      8 1833 1185 (N13W63 L=154)   548     60        0 (1808/1158-1857/1208) -31.0 W59-W65 N11-N14
      9  358 1151 (N08E47 L= 45)   385    161        0 ( 304/1089- 388/1189) -40.6 E52-E43 N04-N10
     10 1749  856 (S07W52 L=144)   509     28        2 (1726/ 831-1772/ 880)  14.0 W50-W54 S08-S05
     11  900  750 (S14E07 L= 84)  1191    280        3 ( 856/ 642- 945/ 876)  78.7 E10-E04 S21-S06
     12 1042  839 (S08W02 L= 93)   280    150        1 (1003/ 819-1089/ 866)  26.6 W00-W04 S09-S06
     13  451  785 (S14E39 L= 53)   266     70        1 ( 425/ 766- 488/ 806)  11.3 E40-E35 S15-S12
     14  182  744 (S19E71 L= 20)   740     70        0 ( 153/ 721- 215/ 760)   9.5 E82-E64 S20-S17
     15 1200  667 (S19W12 L=104)   377    103        1 (1165/ 626-1235/ 709) -45.0 W09-W15 S21-S16
     16 1469  522 (S29W35 L=126)  1685    244        2 (1429/ 434-1521/ 640) -59.0 W31-W42 S35-S20
     17  775  601 (S25E16 L= 75)   525    171        2 ( 691/ 585- 836/ 626) -11.3 E21-E11 S25-S22
     18 1674  495 (S31W57 L=149)   571     88        1 (1632/ 483-1715/ 511)  -9.5 W51-W63 S32-S29
  

    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